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Introduction: The initial accretion of primitive bodies 
(asteroids and TNOs) from freely-floating nebula particles re- 
mains problematic. Here we focus on the asteroids where 
constituent particle (read “chondrule") sizes are observation- 
ally known; similar arguments will hold for TNOs, but the 
constituent particles in those regions will be smaller, or will be 
fluffy aggregates, and are unobserved. Traditional growth-by- 
sticking models encounter a formidable “meter-size barrier" 
[ 1 ] (or even a mm-cm-size barrier [ 2 ]) in turbulent nebulae, 
while nonturbulent nebulae form large asteroids too quickly 
to explain long spreads in formation times, or the dearth of 
melted asteroids [3]. Even if growth by sticking could some- 
how breach the meter size barrier, other obstacles are encoun- 
tered through the l-10km size range [4], 

Another clue regarding planetesimal formation is an ap- 
parent 100 km diameter peak in the pre-depletion, pre-erosion 
mass distribution of asteroids [5]; scenarios leading directly 
from independent nebula particulates to this size, which avoid 
the problematic m-krn size range, could be called “leapfrog" 
scenarios [ 6 - 8 ], The leapfrog scenario we have studied in de- 
tail involves formation of dense clumps of aerodynamically 
selected, typically mm-size particles in turbulence, which can 
under certain conditions shrink inexorably on 100-1000 orbit 
timescales and form 10 - 100 km diameter sandpile planetesi- 
mals. The typical sizes of planetesimals and the rate of their 
formation [7,8] are determined by a statistical model with prop- 
erties inferred from large numerical simulations of turbulence 
[9], Nebula turbulence can be described by its Reynolds num- 
ber Re = ( L/r /) 4 ^ 3 , where L = Ha 1 ^ 2 is the largest eddy 
scale, H is the nebula gas vertical scale height, and a the neb- 
ula turbulent viscosity parameter, and r/ is the Kolmogorov or 
smallest scale in turbulence (typically about 1 km), with eddy 
turnover time t v - In the nebula, Re is far larger than any 
numerical simulation can handle, so some physical model is 
needed to extend the results of numerical simulations to nebula 
conditions. 

Cascade model and multiplier distributions: The spa- 
tial distribution of particle concentration can be captured statis- 
tically by a cascade model [7,9] which predicts the probability 
distribution functions (PDFs) for dense particle clumps; these 
PDFs are essentially volume fractions in the nebula having 
the necessary properties (solids mass and local vorticity) for 
planetesimal formation. A cascade model presumes that, as 
energy flows from large eddies to smaller ones, particles and 
fluid properties are partitioned unequally at each “level" of 
the cascade from “parent" into “daughter" eddies. The proba- 
bility distribution functions (pdfs) of the so-called multipliers 
by which particle and fluid properties are partitioned have 
widths that are given by 1//3; smaller (3 means larger width, 
or a higher probability of strongly asymmetrical partitioning 
which, repeated over a number of levels, leads to clumpier or 
more intermittent local values of particle and fluid properties. 

The wide range of scales between L and r/ is loosely 
called the inertial range; in this range the equations of mo- 
tion are scale-free, and the multiplier pdfs for dissipation of 
turbulent kinetic energy e, which, like our preferred particles, 
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Figure 1: The symbols indicate the multiplier pdfs 

for turbulent dissipation in high-77e atmospheric turbu- 
lence [10], showing scale invariance over 300-300077- The 
red /3-functions (/3=3, 6 , and 48) provide context for (3 
values in figure 2. In [7] we adopted level-independent 
multiplier pdfs with f3 ~ 3 based on our own simula- 
tions. 


is characterized by Kolmogorov eddy timescales, are known 
to be scale-independent in atmospheric flows across hundreds 
to thousands of Kolmogorov lengthscales [10; figure 1], This 
fact motivated us to assume in [7] that the multiplier pdfs for 
particle concentration , which we denote by 4?, were also in- 
dependent of nebula turbulent eddy lengthscale. While this 
assumption was checked and seemingly validated (see figure 
4 of [9]), we have always been concerned about it [7], and dis- 
crepancies between our model predictions and recent results 
of [ 1 1 ] have led us to determine more precisely the range of 
scales over which our assumption is or is not valid. 

Here we report a scale dependence of the widths of the mul- 
tiplier pdfs for both e and 4>, over a range of the largest scales, 
covering a decade or more between the energy-containing or 
integral scale L and the top of a more rigorously scale-free 
inertial range. This scale dependence probably explains the 
differences between [7] and [11], and will affect the formation 
rate and masses of primary planetesimals in the scenario of 
[7,8]. 

Method and results: We obtained a complete dataset 
of particle trajectories in direct simulations covering 6 large 
eddy turnover times (Tff in high Re isotropic turbulence, for 
2 — 6 x 10 s “inertial" particles at each of a range of gas drag 
stopping times ( t s > 0 ), and for fluid tracers ( t B = 0 ) [12, 13]. 
The original simulation was on a cube of 2048?/ on a side, and 
had an integral scale of about 102477. The data are thus similar 
to [9,1 1] but at a considerably higher Re, and showed a clear 
inertial range [12] extending over tens to nearly IOOO 77 based 
on the second order velocity structure function. The particle 
locations were spatially binned into snapshots in time, initially 
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Figure 2: The width of multiplier pdfs for particle con- 
centration is the inverse of the j3 value, plotted here as 
a function of particle Stokes number St = t 3 /t v . Each 
curve represents the situation at a different lengthscale 
(in units of 77). Note that particles with St > 1 but still 
on the order of unity approach /3 ~ 5 — 10 by 6477 ~ L/16 
and generally stay there with the minimum /3 shifting 
slowly to smaller St at smaller scales and deepening 
towards the asymptotic value f3 ~ 3 (dashed line). 

separated by 77,/ 2 to ensure their independence. We then 
determined the multiplier pdfs in the manner described by [9], 
but now over a range of spatial scales between roughly 1 2-5 1 277 
and a range in St = t s /tri = 0.6 — 70, roughly corresponding 
to 1mm - 10cm particle diameter under nebula conditions in 
the asteroid belt region. The dissipation e was known at the 
locations of the tracers, so its multiplier pdfs could also be 
determined as a function of spatial binning scale. The spatial 
distribution of the St = 0 tracers agreed well with a Poisson 
distribution, as they do not undergo preferential concentration. 

Because the number of inertial and tracer particles in this 
study was smaller than used in [9,11], the average number 
of particles per 3277-scale grid cell (for example) was only 
about 2.5. This led to several small-number-statistics sam- 
pling biases that needed to be allowed for. These biases led 
to artificial intermittency; that is, the widths of the apparent 
multiplier pdfs were larger than their true values. This effect 
becomes worse at smaller scales where we typically find fewer 
particles per bin, even if there are more bins. We ran a num- 
ber of Monte Carlo simulations, and found that the bias was 
predictable and iteratively correctable, depending only on the 
number of particles in a bin and the true width of the multiplier 
pdf. For instance, a very broad pdf (small f3) was not biased as 
much by this effect as was a very narrow one. We also realized 
that smaller bins could be sampled on shorter timescales than 
Tl /2, because of the shorter eddy times and faster refresh rates 
associated with smaller lengthscales; this provided improved 
statistics for the smaller bins. 

Figure 2 implies that, at least for the top decade or so of 
eddy lengthscales, multiplier pdfs are slightly narrower (larger 
(3) than we had assumed and particles of St ~ 2 — 10 are 
more concentrated (ie, have broader multiplier pdfs than those 


for St = 1). By contrast, particles of dm-size and larger have 
narrow multiplier pdfs at all but the very largest scales. Our 
preliminary results for t are similar to the particle results, in 
that the multiplier pdf which is so level-independent over a 
wide range of scales in atmospheric turbulence (figure 1) also 
starts out very narrow at the largest scales and only approaches 
j3 ~ 3 at 6477 or so, roughly 1.5 decades below L. We have 
not yet completed the correction of the e results for bias as of 
this writing. 

Thus, the multiplier pdfs for both e and <F paint a consistent 
picture. In the top decade of lengthscales, turbulent stretching 
and vortex tube formation (see [8] and references therein) have 
not yet reached their fully developed state, so the statistics of 
e and $ of particles with Kolmogorov-eddy stopping times do 
not either. Flowever, both e and 4? multipliers approach the 
asymptotic value of /3 ~ 3 at slightly different but comparable 
scales of ~ L/ 30, and we know from atmospheric turbulence 
that this asymptotic value then continues to hold over the en- 
tire inertial range, at least for e. The important implication is 
that multiplier statistics are not characterized by their scale as 
a multiple of rj but as ?l fraction of L. Observed level inde- 
pendence of e to 3OOO77 in atmospheric flows [10] is consistent 
with their much larger L, such that 30007/ < L/16. 

The cascade model for planetesimal formation by turbu- 
lent concentration [10] can now be modified easily to allow 
for the observed level-dependence in the multipliers (figure 
2), and since scales are resolved in the calculations reported 
here to almost the scale where level-independent pdfs are at- 
tained, future uncertainties in calculated IMFs for planetesimal 
formation will be greatly reduced. Moreover, particles that be- 
come most densely concentrated (are most intermittent, or have 
broader multiplier pdfs) in the range of eddy scales near the 
inertial scale have stopping times perhaps 2-10 times larger 
than we previously assumed, consistent with some previous 
work[ll,14]. At smaller scales, St = 1 particles will increas- 
ingly dominate (figure 2). Because nebula properties such as 
turbulent intensity and gas density are not known to factors of 
order unity, these findings do not qualitatively change the out- 
look for the scenario (note the basic timescale revision found 
by [8]), but will lead to quantitative changes in predictions of 
planetesimal IMF and chondrule size distributions. 
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